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Abstract 

The basic principles of classical and semi-classical theories of molecular optical activity are 
discussed. These theories are valid for dilute solutions of optically active organic molecules. It 
is shown that all phenomena known in the classical theory of molecular optical activity can be 
described with the use of one pseudo-scalar which is a uniform function of the incident light 
frequency oj. The relation between optical rotation and circular dichroism is derived from the basic 
Kramers-Kronig relations. In our discussion of the general theory of molecular optical activity 
we introduce the tensor of molecular optical activity. It is shown that to evaluate the optical 
rotation and circular dichroism at arbitrary frequencies one needs to know only nine (3 + 6) 
molecular tensors. The quantum (or semi-classical) theory of molecular optical activity is also 
briefly discussed. We also raise the possibility of measuring the optical rotation and circular 
dichroism at wavelengths which correspond to the vacuum ultraviolet region, i.e. for A < 150 nm. 

PACS number(s): 33.55.+b and 33.20. Ni 



* E-mail address: afrolov@uwo.ca 
t E-mail address: dwardlaw@uwo.ca 



1 



I. INTRODUCTION 



In this study we discuss the classical theory of molecular optical activity. This theory was 
originally developed for solutions of various optically active organic molecules (see, e.g., jlj], 
Q and references therein). Our analysis begins with the classical theory of optical activity 
based on the Maxwell's equations for electromagnetic fields. Any optical active substance 
is described in this theory with the use of a few phenomenological parameters. The main 
goal of the classical theory of optical activity is to derive some useful relations between these 
parameters in various cases. In general, these parameters also depend upon frequencies and 
relations between such parameters take different forms for different frequencies. We also 
consider semi-classical theory of optical activity of molecules which was originally developed 
by Rosenfeld in [3|. In this theory all molecules are considered as quantum systems, while 
radiation is considered classically. This old theory is still widely used, since it produces a very 
good agreement with many experimental results. In particular, the semi-classical theory of 
optical activity can be used at short and very short wavelengths, e.g., for wavelengths which 
correspond to vacuum ultraviolet. On the other hand, it is clear that the complete theory of 
optical activity can be based only on quantum mechanics of molecules and quantum theory 
of radiation. 

This work has the following structure. Basics of the classical theory of molecular optical 
activity in dilute solutions of organic substances can be found in the next Sections. Here 
we introduce the optical rotatory parameter j3. The four Stokes parameters are defined in 
Section III. These parameters are very convenient to describe quasi-monochromatic light. 
The phenomenon of circular dichroism is described in Section IV. It appears that the two 
fundamental a;— dependent functions (optical rotation and circular dichroism) which can be 
defined for an arbitrary optically active solution can be written in the form of one complex 
function. The well known Kramers- Kronig relation between the real and imaginary parts of 
this functions must always be obeyed. For limited intervals of frequencies this produces a 
very useful relation which allows to determine, e.g., the circular dichroism by using the known 
expressions for optical activity. Tensor of molecular optical activity is explicitly defined in 
Section V. The formulas obtained in this Section are very useful in applications, since they 
allow to express the optical activity by using only the two basic molecular properties (the 
electric dipole and magnetic dipole momenta). Concluding remarks can be found in the last 
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Section. 



II. CLASSICAL THEORY OF MOLECULAR OPTICAL ACTIVITY 

Let us briefly discuss the classical theory of optical activity. In classical theory the op- 
tical activity always denotes the ability of the material under study to rotate the plane of 
polarization of the left- and right-circularly polarized light. Currently, the study of opti- 
cal activity also includes optical rotation at different wavelengths, circular dichroism, and 
differential scattering of left- and right-circularly polarized light. All these phenomena are 
manifestations of natural optical activity which is a characteristic of chiral molecules (in 
contrast with achiral or non-active molecules). Note that there are also various phenomena 
which correspond to so-called induced optical activitivies. In such cases the achiral molecules 
can show some optical activity, if, e.g., they are placed in a relatively strong electric and/or 
magnetic fields. In this study we restrict ourselves to the analysis of the natural optical 
activity only. 

In general, the optical activity is uniformly related to the spatial dispersion, i.e. to the 
non-local relation between the electric induction D and electric field E. For the Cartesian 
components of these vectors we can write {4] 

/*oo P 

AO, t) = Ei(r, t) + J drj d^F^r, r, rjE^t - r, n) (1) 

where Fik(r, r, r x ) = F^ir, r l5 r) is the kernel of integral operator. For monochromatic field 
components E(r, t) = E(r) exp(— iu)t) and Eq.flT]) takes the form 

A(r) = #i(r) + J (pTjytoT^Ejfa) (2) 

This equation with the kernel fij{oj; r, ri) expresses a non-local relation between D and E 
which is also called spatial dispersion. In general, the kernel fij(u;r, ri) in Eq.Q rapidly 
decreases with interatomic distances. In many cases such a kernel is very small already at 
distances ~ 3 a, where a designates some average (or effective) atomic dimension. Briefly, 
the relation, Eq.(j2J), is written in the form 

A(r) = J d 3 ri ■ e^u- r, r^Efa) (3) 

where e^- is some non-local operator (tensor) which also depends upon the frequency u (see 
below). The 3x3 tensor e = is the dielectric tensor (also called electric permittivity). 
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From the transparency of the media it follows that all nine matrix elements are real and 
all three eigenvalues of this tensor are positive. 

In all studies of optical activity in organic compounds only infinite, homogeneous media 
are considered, a convention adopted in this work. Unless otherwise specified, the absorbtion 
of radiation is assumed to be absent at all frequencies considered below. In such cases the 
kernel in Eq.([21) depends only on the difference R = r — ri. The functions D and E 
in infinite, homogeneous media can be expanded in a Fourier integral with the respect to 
Cartesian coordinates as well as time. Finally, this allows one to obtain the following relation 
between the corresponding Cartesian components of the vectors D and E 

A(k) = ey(a;; k)%(k) = [<% + jT ' dr J d 3 R^(r, R) exp[i(wr - k • R)]] E^k) (4) 

In other words, the dielectric tensor eik{oo;k) (also called electric permittivity) takes the 
form 



/*oo P 

[u, k) = 5ij + J dr J d 3 Rfij{T, R) exp[i(ujr - k ■ R)] (5) 



As follows from this formula the dielectric tensor is a function of the field frequency u and 
wave vector k. In general, the dependence of the dielectric tensor e^- on u is called dispersion, 
while the analogous dependence upon the wave vector k represents the spatial dispersion. 
The spatial dispersion of Cjj(k) is responsible for optical activity (see below). 

In solutions of organic substances the optical activity corresponds to the case of weak 
spatial dispersion, i.e. k —\ k | is small. In such cases the tensor £y(u;,k) can be expanded 
in powers of the wave vector k, e.g., 

6ij(u, k) = e\f(uj) + Jiji(uj)ki + Piji m (oj)kih m + aiji mn (u;)kik m k n + ... (6) 

Such an expansion is valid, if the first term in Eq.([51), i.e. e^\u), has no zeros in a given 
range of frequencies u. Since in this study we restrict ourselves to the consideration of 
transparent (or slightly absorbing) solutions only, then we can neglect the imaginary part 
of dielectric tensor e\f(u). 

If these two conditions are obeyed, then for small k only the first few terms in such an 
expansion are important. Let us restrict to the two lowest order terms (the second of which 
is responsible for optical activity), i.e. we can write 

ey(w, k) = e\f(u) + Jiji(uj)ki = e^\u) + i-~/ m {u)ni (7) 

4 



where n = and j^ni is an antisymmetric tensor of the second rank (upon indexes i and 
j). For tensor 7^ the antisymmetry means 7^7 = —Jju- If absorbtion of radiation is absent, 
then the tensor 7^7 is real, i.e. 7*7 = 7^7. These two conditions mean that the jijirii tensor 
can be re-written into another form jijini = ^e^igi, where eiji is the complete antisymmetric 
tensor, while g\ is the I— th component of the axial giration vector g. In tensor algebra this 
relation is called the duality relation. In general, the giration vector g is a function of the 
unit wave vector n, i.e. (g)j = guni, where gu is the pseudotensor of the second rank. In 
isotropic media gu = 5uf and such a pseudotensor is reduced to a single pseudoscalar /. In 
this case the 7^7 tensor is essentially reduced to the complete antisymmetric tensor e^. 
In fact, for the tensor 7^7 one finds 7^7(0;) = ^eijif(u), and therefore, 

D = e (0) E + i/(w)(Exn) (8) 

Note again that this equation can be applied only in those cases when the absolute value of 
f(co) is much smaller than the minimal eigenvalue of the tensor e^°\u). 



As is well known (see, e.g., [4]) in an arbitrary dielectric media we always have D ■ n = 0. 
In this case from Eq.flH]) one also finds that E • n = 0. For a monochromatic wave we can 
write the Maxwell equations 

= (k x E) and = (k x H) (9) 

It follows from here that k _L D _L H and also that E 1 H. In three-dimensional space 
this means that the three vectors E, D and k are coplanar. This simplifies drastically the 
following analysis of optical activity in isotropic media. 

Consider now the energy transfer. In general, the direction of the energy flux is given by 
Poynting vector S = ^(E x H). Now by using the unit vector n defined above (n = -k) 
we can write for the Poynting vector 

S = ^-(ExH) = ^K-(E.n)E] (10) 

The total energy flux through an element dS of surface orthogonal to n is 

dW = 7— [-E 2 — (E ■ n) 2 ]dQ = -^-E 2 sin 3 GdQd<5> (11) 

47T 47T 

where is the angle between the vector E and outer normal to this surface element dS, i.e. 
n. Also, it follows from the two equations of Eq.()9]) that D = n 2 E — n(n ■ E). On the other 



hand the basic relation between vectors D and E is D = eE, where e is the dielectric tensor. 
From here one finds the following equation written in Cartesian components 

(n 2 5ij - mrij - eij)Ej = (12) 

where are the components of dielectric tensor. 

Formally, this equation coincides with the corresponding eigenvalue equation for the 
dielectric tensor e^. However, the eigenvalues of this tensor are the functions of three 
spatial directions. By using some unitary transformation, one can reduce Eq.( TT2|) to the 
principal axes of the tensor ey which are also called the principal dielectric axes. In fact, 
there are some advantages to writing Eq. fll2j) in the principal dielectric axes. In this case 
it exactly coincides with Fresnel's equation which is the main equation of crystal optics. 
In general, Eq. (ll2j) determines the wave- vector surface in the n x ,n y ,n z coordinates. Such 
surfaces depend upon three constant coefficients e x , e y , e z (eigenvalues of the dielectric tensor 
e tj ). 

For homogeneous solutions the Fresnel's equation simplifies significantly, since is these 
systems e x = e y = e z . We want to consider such a transition in the two following steps. 
First, consider the case of two different eigenvalues e x = e y = e±_ and e z = eii (these values 
of parameters correspond to uniaxial crystals). In this case the Fresnel's equation can be 
factorized to the form 

{n 2 - eA_)[e\\n 2 z + e ± (n 2 x + n 2 y ) - ej_e||] = (13) 

where n = (n x , n y , n z ) is the direction of the light propagation. In other words, an equation 
of the fourth order (upon n) is reduced to the product of the two quadratic equations 

n 2 = n 2 x + n 2 y + n\ = e± (14) 
^ + ^4 = 1 (15) 

where the first equation is the equation of a sphere, while the second equation determines 
an ellipsoid. The sphere represents the propagation of ordinary waves. Such waves have the 
same refractive index n = y/e±. The second equation represents the so-called extraordinary 
waves which are directly related with the optical activity. Let us consider the extraordinary 
waves in homogeneous solutions, or in crystals of a cubic system. These two cases can be 
obtained as the limit of Eq. f|T5|) when eii — > e±. In reality, it can be written in the two 



different forms eii = e± ± A, where the positive parameter A — > 0. In such cases, Eq.flT 
does not change, while the second equation takes the form 

n 2 n 2 + n 2 n 2 cos 2 6 n 2 sin 2 9 , , 

e± e±±A e± e ± ± A v ' 

where n x = n ■ sin 9 cos (p,n y = n ■ sin 6 sin (p,n z = n ■ cos 9, where 9 is the angle between the 
optical axis and vector n. 

In homogeneous solutions the orientation of chiral molecules is random, i.e. we have to 
replace the factors cos 2 9 and sin 2 9 in the last equation by their mean values ~. This gives 

^ + drA = l (17) 

or 

2 _ e ± + e ± ± A 2(e ± ±]A) = 2 

^ e ± (e ± ±A) ~ (e ± ± |A) 2 (e ± ± |A) 1 J 

From here one finds that n 2 = e± ± |A, or in other words, we have two different refractive 
indices n\ = e_|_ + |A and n\ = e± — |A. This means that two different refracted wave are 
formed and, formally, we have to consider the double refraction or birefringence. However, 
the parameter A is small (in fact, very small) in comparison with n 2 . Therefore, the overall 
scale of such a birefringence is ~ A. 

Result can be obtained in a slightly different way with the use of some microscopic 
identities. Indeed, let us note that for homogeneous solutions of chiral molecules = 
e ■ 5ik + i~f(u))en~ini, where e^i is the complete antisymmetric tensor. In this case we do 
not need to use the complete Fresnel's equation to produce the same answer as above. The 
chiral activity can be described with the use of only one numerical parameter f(oo) which 
is pseudoscalar. In fact, such a parameter can be introduced in a slightly different way. 
Indeed, the Maxwell equations in the case of homogeneous solutions take the form 

OH <9E 
B = eE-g— , B = uH + g— (19) 

where g is a some constant, which depends upon u. By taking into account polarization of 
media by the electric and magnetic field we can write, e.g., for the D vector 

D = (1 + ATcNa)E — 4nN——- (20) 

c at 

where a is the static polarizability, while /3 is the so-called optical rotatory parameter, or 
optical rotation, for short. As follows from Eq. (l2"U|) the parameter j3 is also a pseudoscalar. 



The parameter /3 determines the optical rotation, i.e. the rotation of the plane of left- and 
right-circularly polarized light when it passes through an optically active medium. Also, in 
this equation iV is the number of chiral molecules per unit volume. It follows from the last 
two equations that e = 1 + AirNa and g = AtcN^. The relation between factor g and the 
indices of refraction for circularly polarized light can also be obtained from Eq. (l20l) 



Xl = y/e - 2nug and XR = \/t + ^ug (21) 

Now it is easy to find the overall rotation (5) when the light propagates the distance z in 
some chiral media 

irz , . Att 2 z A 9 o 167T 3 i> 2 Nz n 
S = -j(xr - Xl) = —g = 4tt 2 zA ■ g = (22) 

Note that the optical rotatory parameter f3 (as well as a) which follows from Eq. (1201) can 
rigorously be determined only with the use of the quantum mechanics of molecules. This 
will be our goal in the third Section. 

III. STOKES PARAMETERS 

As follows from its definition any monochromatic wave has a certain polarization. How- 
ever, in actual optical experiments it is almost impossible to create a beam of pure monochro- 
matic waves, and usually we have to operate with real light which contains frequencies dis- 
tributed in a small interval Au around the main frequency u. The means that the real light 
is, in fact, a mixture of light quanta with different polarizations. An arbitrary property of 
such a beam of light, e.g., the electric field E, in the real light depends upon time. If the 
frequency distribution Au around u is narrow, then the E(t) function can be represented in 
the form E(t) = E (£)exp(— tut), where E (£) is a slowly varying function of time t which 
determines the polarization of the actual light. From the last formula one can expect that 
such a polarization will be slowly changing in time, i.e. we are dealing with the partially 
polarized light ^]. 

In regular experiments we cannot observe the polarization properties of electromagnetic 
waves directly. Instead, one measures the intensity distribution of light when it passes 
through various physical bodies. This means we are dealing with quadratic functions of 
the field. In other words, in actual experiments we are measuring the components of the 
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following tensor J a p = E 0o ,Eqo, where E 0a and Eqq are the Cartesian components of the 
slow varying E (t) vector. The line over the product of the two complex vectors mean the 
value averaged in time. If all vectors are represented in the form E(t) = E (t)exp(— toot), 
then the Eo a EQg product is the only value which has non-zero time-average. Other similar 
combinations, i.e., E^E^a and E 0a E 0/3 , contain rapidly oscillating factors exp(—2iut) which 
gives zero upon time averaging. 

Since in any plane wave one finds El n, where n is the direction of wave propagation, 
then the J a p tensor has only four components. Moreover, the J a p tensor also contains the 



total intensity of the wave J = Y, a Jaa — E • E*. This value has nothing to do with with 
polarization of the wave and can be excluded by introducing the tensor p a $ = - s j-. The 
tensor p a p has the unit trace and it is called the polarization tensor. It can be shown that 
the polarization tensor is hermitian, i.e. p*^ = pp a . Now, we can introduce the degree of 
polarization P which is defined as 



P = yjl - Adet(p aP ) = yjl - Ap n p 22 + 4 | p 12 | 2 (23) 

where det(pp a ) is the determinant of the 2x2 matrix pp a . 

An arbitrary hermitian 2x2 matrix can be represented in the following form 

1 1 i 

Pafi = ^(pa/3 + P/3a) + ^(/V ~ P/3a) = Sa/3 ~ ^e a/3 A (24) 

where S a/ 3 is the real symmetric 2x2 tensor. The analogous non-symmetric 2x2 ten- 
sor in two-dimensional space is reduced to the unit antisymmetric tensor ei 2 = —e 2 \ and 
pseudoscalar A. The pseudoscalar A is called the degree of circular polarization. It is 
bounded between -1 and +1. The case of a light wave with linear polarization corresponds 
to the A = value. The waves with circular polarization correspond to the values A = +1 
(right-circular) polarization and A = — 1 (left-circular polarization). 

An alternative expansion for an arbitrary hermitian 2x2 matrix is performed with the 



use of three Pauli matrices <7j (i = x, y, z) (see, e.g., jsj) and one unit matrix /. It is written 
in the form 

1 ( 1 + £:. (.1 - <2 

2 



PafS = ~ {I + £l<T w + fay + J = 9 ( 25 ) 



The parameters ^1,6 and £3 which appear in this formula are the so-called Stokes parame- 
ters. In general, any intensity measurement may be written as a linear combination of these 
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three parameters and one additional Stokes parameter £o which is the total intensity of the 
scattered light. The determinant of the pp a tensor is 



while the degree of polarization is P = y£f + £2 + £!• The Stokes parameters £1 and £3 
determine the degree of the linear polarization, while the parameter £2 shows the degree 
of circular polarization. Note that the parameter £ 2 coincides with the pseudoscalar A 
introduced above. From three Stokes parameters one can construct the two scalars (£2 and 



£1 + £3) which are invariants under Lorentz transformations. The three Stokes parameters 



IV. CIRCULAR DICHROISM 

In all formulas above we have neglected the absorbtion of light during its propagation in 
the dense media. In actual cases the absorbtion of light always occurs. At some frequencies, 
e.g., in the vacuum ultraviolet region, it plays a very important role and cannot be ignored 
even in the first approximation. In reality the situation is even more complicated, since 
light waves with different circular polarization are absorbed differently by the media. This 
is called circular dichroism (CD). Such a differential absorbtion of light with left- and right- 
circular polarizations can directly (and substantially) affect the observed optical activity. 
It appears that optical rotation and differential absorbtion of light with different circular 
polarizations can be considered as the two manifestations of one phenomenon. 

In general, a detailed study of circular dichroism at different frequencies allows one to 
develop a new approach for analysis of organic substances. In this Section we want to discuss 
the modification which is required in all formulas presented above. The absorbtion of light 
is described by introducing an imaginary part into the permittivity tensor €ij(u,r), or in 
other words, by considering static polarizability a as a complex value. However, we are 
not interested here in the total light absorbtion. Our interest is related to a very specific 
difference between absorbtion of light with left- and right-circular polarizations. It is clear 
that a complex static polarizability a, Eq. (j20JL cannot describe such differences. As follows 
from Eq. (1201) this goal can be achieved by considering the optical rotatory parameter f3 as a 
complex value. 




(26) 





also have a number of other advantages in actual applications [7|. 
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In these cases the parameter is represented as the sum of its real and imaginary parts, 
i.e. /3 = fti + ifa, where /3% and are two functions of the frequency u. These two functions, 
however, are not completely independent, since there are two additional connections between 
them which follow from the Kramers- Kronig relations. This follows from the fact that (3(uj) 



is the response function [4|, J9] which is an analytical function in the upper half uj plane (for 
now, we consider the frequency to as a complex variable) (see, e.g., [9|). This allows us to 
use Cauchy's theorem for (3(u): 

The contour C can be chosen to consist of the real frequency axis uj and a great semicircle 
at infinity in the upper half plane. The function f3(u) vanishes rapidly at infinity, i.e. there 
is no contribution to the integral from the great semicircle. Finally, Cauchy's integral is 
written in the form 

w+~^w (28) 

where z now is any point in the upper uj— half plane and the integral is taken over the real 
axis. In fact, we want to place the point z at the real axis. This can be done by approaching 
the real axis from above, i.e. by representing the complex variable z in the form z = uj + is. 
This gives 

2m J-oo Uj' — UJ — IE ,-, 



The denominator in the last formula can be written in the form 

1 



9] 



- = p( — - — ) + ni6(u' - u) (30) 

where the symbol P means the principal part, while 8{x) designates the Dirac delta-function. 
Now Eq.( l29|) takes the form 

Lp f^^W (31) 

Til J-oo UJ — UJ 

By separating here the real and imaginary parts one finds 

^ n , s l n r +OQ ImB(uj')duj' 

Re/3(uj) = -P / (32) 

7T J-oo UJ — UJ 

i mm = - i - P r Re ^' w 

7T J-oo UJ' — UJ 

This is the most general Kramers-Kronig relations written for the optical rotatory parameter 
(3. In general, it can be shown that the Rej3{uj) is an even function in u, while Imf5{uj) is 
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odd. This allows one to transform the last two integrals in Eq. (l32j) to the integrals taken 
over positive frequencies only, i.e. 



These formulas can be used in actual applications which include the optical rotatory pa- 
rameter 0. By using Eq. (122|) we can re-write these formulas for the corresponding angles 

5 = + in 



The importance of the linear frequencies v follows from the fact that these values are usually 
used in actual experiments. In general, Eq. flMJ) represents the explicit relation between the 
actual optical rotation (angle 9) and circular dichroism (angle k). As follows from Eq. (|34p 
the optical rotation known for all frequencies allows one to determine the circular dichroism 
at each frequency In reality, however, one finds a number of restrictions which exist 
in the solution of this problem. Most of such restrictions follow from the fact that optical 
rotations in the VUV region (A < 150 nm) are not known even approximately. On the 
other hand, it is clear that for each molecule the VUV area of wavelengths contains many 
resonance lines which are crucially important to describe the absorbtion of radiation. If 
we ignore the VUV region of wavelengths, then we can restore the circular dichroism at all 
frequencies only approximately (in reality, very approximately). For some limited areas of 
wavelengths, however, such a reconstruction can be quite accurate and complete. Usually, 
these areas of wavelengths are located far from the VUV region. 

Note also, that the experimental knowledge of the 9{v) and k(v) values for large number 
of different frequencies z^, . . . , v n is used to detect uniformly the corresponding organic 
substance. Formally, such an identification allows one to solve many problems of quantative 
and qualitative analysis of the mixtures of chiral organic substances. 




(33) 




(34) 



where we have also introduced the linear frequency v = ^ 



(oo is called the circular frequency) . 
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V. TENSOR OF MOLECULAR OPTICAL ACTIVITY. ROTATION POWER. 



In the middle of 1930's Placzek shown 10] that a significant number of effects related to 



the interaction between atom(s) and electromagnetic field can be described with the use of 
only one tensor, later known as the tensor of light scattering. In particular, the differential 
scattering cross-section of light by an atom (or any other electron containing system) can 
be written in the form 

da = "(" + ^) 3 | (C lk ) 2l {e'-ye k | 2 do' (35) 

it c 

where (6^)21 is the 3x3 tensor of light scattering, while e« and are the polarization 
vectors of the incident and final photons. The integration in Eq. (l35j) is performed over 
the angular variables of the final photon which is designated by the superscript '. Here 
and everywhere below we shall assume that the angular volume element do' has the form 
do' = sin9' d9' d(f)' . The explicit expression for the light scattering tensor (6^)21 is H] 

(C-A _ y^l (dj)2n(dk)nl j (dk)2n(dj) n i - 

„ \-u n i - u — ^0 u n i + to' — J 
where u' = u> + Uu, while di and dk are the corresponding components of the vector of 
the dipole moment d. Note that the differential cross-section do, Eq. fl35|) . corresponds to 
the lowest order approximation upon the fine structure constant a ~ ^ and contains only 



the electric dipole-dipole interaction. The method developed by Placzek [10] was based on 



eqalier studies by Kramers and Heisenberg 12| and Dirac [13 1. 



The Placzek approach for atoms suggests attempting to derive an analogous method 
for molecules which would describe their optical activity. In this Section this problem is 
considered in detail and it is shown that in the lowest order approximation can be described 
by the tensor (6*^)21 of light scattering and by the four (or two in some cases) new tensors. 
These tensors are called the tensors of (molecular) optical activity. Note that in many 
actual cases the four/two tensors of optical activity are reduced to one tensor only. To 
produce the closed analytical expressions for these tensors below we shall assume that the 
electromagnetic field is represented as a combination of plane waves. Each of these plane 
waves has its own frequency 00 and polarization which is represented by the vector e. The 



wave functions of the incident and final photons can be taken in the form (see, e.g., 11| ) 

A ew = \ — exp(—iojt + zk • r)e , A e 'aj' = \ —exp(—icu't + ik' ■ r')e' , (37) 
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where u and u' are the corresponding frequencies, while vectors e and e' represent the 

polarization of the incident and final photons, respectively. Below, we shall consider the 

plane waves in the transverse (or radiation) gauge, where divA = 0. In this gauge one finds 

k ■ e = and k' ■ e' = 0. Note that in calculations for the final photon we need to use 

the wave function which is conjugate to its wave function, i.e. A*, w ,. As follows from these 

equations the electric E and magnetic H fields are 

d 

E ew = ——A eul = —iv 27tujeexp(—iujt + ik ■ r) (38) 

[27T 

H eu] = curlA eul = id — (k x e)exp(—iujt + ik ■ r) 
y to 

By introducing the unit vector n = - we can re-write the last equation in the form 



H ew = iv2tujj{t\. x e)exp(—iu>t + ik ■ r) (39) 

Analogous expressions can be obtained for the E e '^/ and H e y fields which are related with 
the A*, u , wave function. 

From these equations one finds the following expressions for the electric dipole and mag- 
netic dipole interactions. In fact, for each of the (e, o;)-components of the E and H vectors 
we have 

V£, = -d ■ E ew = iV2nuj(d ■ e)exp(-iut + ik ■ r) (40) 

and 



V™ = — m • H ew = -iV2~7Wj[m ■ (n x e)]exp(— uvt + ik ■ r) (41) 

where d and m are the vectors of the electric and magnetic dipole moments, respectively. 
In the lowest order approximation the one-photon matrix elements of the V e and V m in- 
teractions equal zero identically. The first non-zero contribution can be found only in the 
second order of perturbation theory. In the second order approximation the matrix element 



V21 for the transition between states 1 and 2 is written in the following form 11] 

n °1 °n °1 °ri 

where the notation 8 designates the total energy of the system ('molecule + photons'), i.e. 
in the case considered here we have = E n and £^ = E n + u + u'. The matrix elements 
V a b represent absorbtion of the photon with the wave vector k. Analogously, the matrix 
elements V ah represent emission of the photon with the wave vector k'. In the general case, 
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in Eq. (|42|) the V = V 2 i interaction is represented in the form V = V e + V m + V qe + V qm + . . . 
and V = iV e )' + iy m )' + (V qe )' + (V qrn y + where V e , V m , V qe are the electric dipole, 
magnetic dipole and electric quadruple interactions, respectively. Keeping only lowest order 
terms in the expansion of V in terms of the fine-structure constant a ~ we can write 
1/ « 1/ e + V m and V w (F 6 )' + (V m )' . In this case one finds from Eq.(g2I) 



t / V ;2n"nl | V v ;^!nV K Vnl , V v ;2n r nl , V v Jnl , V v V2n " nl //i o\ 



P. —PII I I C c?/ + 



+ 



— £11 \ ' c _ pi p pll 

By neglecting here by all terms ~ ymym an( j g^hgr terms of higher orders in the fine 
structure constant a, we obtain the following formula for the differential cross-section of 
light scattering da 

do - =1 V21 T 4 2 = rf o ee + da em (44) 

where da ee is the part of the total cross-section which can be reduced to the expression given 
above (see Eq. fl35l) ). This part of the cross-section is not related with the optical activity. 
The second term in the right-hand side of Eq. (|44|) is significantly smaller, in the general 
case, than the first term, i.e. da em <C da ee . However, the second term in Eq.f j44l) is a great 
interest, since it represents all lowest order effects which are determined by the molecular 
optical activity. 

As follows from Eq. ( 143)) in order to determine the part of the total cross-section respon- 
sible for molecular optical activity in the lowest order approximation we need to obtain the 
explicit formulas for the matrix elements of the V e V e ,V e V m , and V m V e products. The 
arising expressions are extremely complicated, since each of the V e and/or V m interactions 
contains an infinite number of and V™ components. In the V e V m and/or V m V e prod- 
ucts one finds an infinite number of cross-terms which explicitly depend upon coordinates. 
These terms cannot be computed without a complete and accurate knowledge of the molec- 
ular electron density p e ( r )- 

However, we can introduce an approximation that the wavelengths A of the incident 
and final photons are significantly larger than typical linear sizes of molecule a (our light 
scatterer). In this case we have k-r<|k||r|^C^^0. In this approximation one finds 



from Eqs.dH and QSJ 

E ew = — %\J 2nujeexp(—iujt) (45) 
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and 

H ew = t\ / 2nu(n x e)exp(-iut) . (46) 

Therefore, we can write 



= iV2mj(d ■ e)exp(-iujt) = tV2nu(d UJ • e) (47) 

and 

V™ = -t\ / 2nuj[m ■ (n x e)]exp(—iut) = -iV27ru~j[m 0J ■ (n x e)] (48) 

where d u and m w are the corresponding Fourier-components of the dipole and magnetic 
moments of the molecule. Note that with the identity m ■ (n x e) = e ■ (m x n), the formula 
for magnetic interaction can also be written in another form 

VZ = -iV2^[(m„ x n) ■ e] (49) 

which is similar to the formula for in which the dipole moment d is replaced by the 
vector-product m w x n. 

Now, we can write the lowest order term upon the magnetic interaction in the differential 
cross-section da of the light scattering 

do _ _ I (d 2n ■ e')(d n i ■ e) (d 2ra ■ e)(d ra i ' e ') II ( d 2» ; gOK^ni x n) • e] 

I „ UJ n \ — UJ — lO U n \ + U)' — lO iO n i — id — lO 



[(m*,„ x n) • e'](d„i • e) (d 2n • e)[(m T * tl x n) • e'] [(m 2n x n) • e](d nl • e') u(u') 



+ _ — L + 



A3 



-do' 



U) n i — UJ — lO U n \ + UJ 1 — lO U n \ + UJ 1 — zO 

where the notation e' designates the vector (e')*. This notation is also used in the two 
following equations. This equation can be re-written as 

do _ = I (d 2n ■ e')(d nl ■ e) (d 2ra ■ e)(d ra i - e ') II ( d 2» ; e')[m nl ■ (n x e)] 

I „ U nX — UJ — lO U n i + U' - lO ' ' „ UJ nl - UJ - lO 



[m* 2n ■ (n x e')](dm • e) (d 2n • e)[m* nl • (n x e')] , [m 2n • (n x e)](d ni • e') u(u 



+ — /L , - + 



/\3 



-do' 



U n i — UJ — lO UJ n i + Uj' — lO UJ n i + UJ' — lO 

In these equations and below the notation m* stands for the vector which is a complex 
conjugate vector to the vector of magnetic dipole moment m. In quantum mechanics (in 
the coordinate representation) we always have d* = d, but m* ^ m. Note that the vector n 
in these equations corresponds to the direction of the scattered light. Formally, this vector 
can be oriented in an arbitrary spatial direction, but in almost all modern experiments 
on optical activity in homogeneous solutions the direction of the scattered light always 
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coincides with the direction of the incident light. This means that our differential cross- 
section must be multiplied by a delta-function 5(n in — n) and integrated over the angular 
variables d = (9' , 4>') of the unit vector n = (cos9' cos<p' , cosO'sincf)', sinO') which represents 
the direction of the final photon. The unit vector n in describes the direction of the incident 
photon. This produces the following expression for the cross-section a 

Attu(uj + u 12 ) 3 I ^ (d 2 „ ■ e')(d„i ■ e) (d 2n ■ e)(d nl ■ e') n ^ (d 2n ■ e')[m nl ■ (n in x e)] 



a 



h 2 c 4 I „ w„i — u; — zO oj n \ + u' — lO 1 1 „ u; n i — ui — tO 

[m* 2n ■ (n in x e')](d n i ■ e) + (d 2n ■ e)[m* nl ■ (n in x e')] + [m 2w ■ (n in x e)](d ra i ■ e'jh ^ 



W n l — — lO U} n i + Lj' — lO U n \ + OJ 1 — lO 

where u' = u+Ui 2 and unit- vector n in designates the direction of propagation of the incident 
photon. 

The expression, Eq. f !52j) . can be cast in the following form 

4ttuj(u + U12) 3 



a 



(C ifc ) 2 i(e^)*e fc • (S ik ) 2 i(e )*(n in x e) k + (T ik ) 21 (n in x (e )*);e fe (53) 
+ (Uik) 21(e) i(n in x (e')*) fe + (V ik ) 2 i(n in x e')j(e' fe ) 



where (S ik ) 2 i, (T ik ) 21 , (U ik ) 2 i and (V ik ) 2 i are 3x3 tensors, while the dipole-dipole tensor 
(C ik ) 2 i is defined above in Eq.f l36|) . Here we assume that, in the general case, the vectors 
e' and e which represent the polarization of light are complex. Each of these tensors is 
represented as a sum of its irreducible components, e.g., Si k = S°5i k + S^ k + Sf k , where 

S° = ' ^ik = 7^(Si k + S k i) — S°5i k , S ik = —(Si k ~~ S k i) (54) 

Note also that S°,T°,U° and V° are called the scalar parts of the S,T,U and V ten- 
sors, respectively. The components with the superscripts s and/or a (e.g., T s ,T a ) 
are the symmetric and antisymmetric parts of the tensor. All components of the 
S°, T°, U°, V°, S? k , Tf k , U? k , Vf k , S? k , Tg, U? k and V& tensors contain the products of the cor- 
responding components of the d and m vectors, which are the vectors of the electric dipole 
moment and magnetic dipole moment, respectively. The vector-operators which represent 
the electric and magnetic dipole momenta are assumed to be self-conjugate. Furthermore, as 
mentioned above in the coordinate representation the vector d is a real vector (i.e. d* = d), 
while the vector m is a complex vector (i.e. m* ^ m). For instance, the explicit expressions 
for the S°,T°, U° and V° tensors (they are called the scalar-components) are 

(Q0\ 1 V (^)2«(™i)«l /^ON 1 ( m i)2n(di)nl 

<J n U nl — U <J n W n \ — U 
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mO\ 1 (^i)2n(^j )nl / t/ qs _ 1 (mi)2n(4)nl 

respectively. Analogous formulas for the symmetric and antisymmetric parts of the S, T, U 
and V tensors are significantly more complicated. These formulas and the physical meaning 
of all irreducible components of these S, T, U and V tensors will be discussed elsewhere. 

Thus, we have shown that all phenomena related to the optical activity can completely 
be described with the use of only four tensors: S 2 i, T 2 i, U21 and V 2 i- The fifth ten- 
sor C21 (the tensor of electric-dipole light scattering) is included in the formula for the 
cross-section as an amplification factor. These five tensors have fifteen irreducible tensor- 
components C°, C? k , C? k , S°, T°, U°, V°, T* k , Sf k , U? k and V& The first three 

tensors C°,C- k ,Cf k here have nothing to do with the optical activity itself. Instead they 
determine the amplification factor which also appears to be u— dependent. The optical 
activity is described by the twelve tensors (S°, T°, U°, V°, S? k , T? k , U? k , Sf k , Tg, Uf k and 
V ik ). In many real applications, however, the total number of independent tensors can be 
reduced. For instance, if the 1- and 2-states are identical and u 2 i = (Rayleigh scattering), 
then to describe optical activity one needs only two tensors (not four!) with six irreducible 
components. This case corresponds to the regular optical activity (optical rotation) mea- 
sured in modern experiments with dilute solutions of organic molecules. Furthermore, if the 
polarization vectors are always chosen as real (not complex), then to describe the optical 
activity one needs only one 3x3 tensor with three irreducible components. However, the 
explicit u— dependence of such a tensor will be quite complicated. All such cases will be 
considered in our next study. 

The intensity of the scattered light I' is uniformly related to the intensity of the incident 
light I by the relation 

r=r^W (56) 



u ■ 

As follows from the formula for the cross-section a, Eq. (l52|) . in any optically active solution 
the intensity of the (scattered) light will always be rotated during its propagation along the 
direction n in . The factor (77) f in the last formula can be considered as the rotation power. 
As follows from the last formula the uniform combination of the twelve tensors mentioned 
above multiplied by the amplification factor, Eq. (l36l) . allows one to determine the so-called 
rotation power of the given optically active solution for the initially polarized light. Note 
that only our approach produces the correct and complete formula for the ui— dependence 
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of the rotation power. 



VI. QUANTUM THEORY OF MOLECULAR OPTICAL ACTIVITY 

As can be seen above, the physical origin of the relations between different parameters 
used in classical theory of optical activity remains unknown. The corresponding analyti- 
cal expressions, numerical values and all possible relations between such 'phenomenological' 
parameters can be found only with the use of modern quantum theory based on Quantum 



opment of quantum origin 
3]. A simplified and com- 



14|. Below, 



Electrodynamics. The first and very important step in the deve 
of optical activity was made by Rosenfeld almost 80 years ago 
plete version of the approach developed by Rosenfeld can be found, e.g., in 
we shall follow the same general direction. Our main goal in this Section is to obtain the 
relation between the optical rotatory parameter j3 from Eq. (|22|) and properties of an isolated 
molecule. For an isolated molecule the optical rotatory parameter /3 can only be a function 
of the molecular 2 £ -pole moments. In reality, however, only a few moments with small £ 
(£ = 1,2) contribute noticeably. As we mentioned above the parameter /3 is a pseudoscalar. 
Therefore, the first (largest) term in the expansion of j3 in terms of 2^-pole molecular mo- 
ments is proportional to the scalar product of the dipole vector and the pseudovector of the 
magnetic moment d • m. The second term must be proportional to the product d • Q ■ m, 
where Q is the second order tensor of the electric quadropole moment. 

Below, we shall assume that all molecular wave functions (for the ground and excited 
states) are known (or can be determined) to very good accuracy. In this case, by using 
Rosenfeld's formula one can calculate the optical rotatory parameter (3 (in some studies it 
is also called the chiral response parameter) 



Im 



(a | d | b) (b | m | a) 



p=-a*Y. 1 —* _ , 2 J (57) 



67r ^ T V ab - V 2 

where the summation is taken over all intermediate states. In this equation we use the linear 
frequencies v instead of circular frequencies u, where u) = 2iru. The notation Im designates 
the imaginary part of the terms written in brackets. Symbols a and b stand for the quantum 
(molecular) states, while | a) and | b) mean the corresponding wave functions. Rosenfeld's 
formula is based on an assumption that all molecular states (ground and excited) have zero 
widths. In other words, these states are stable, i.e. the decay time is infinite. In general, 
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this is not a very realistic assumption and we need to introduce finite line widths, ej?;., 
lab(v) = 4?r 6 ^ U 3eahc? a ~ Ub ^ ( m Si-units and in the lowest order dipole approximation 15]). 
The Rosenfeld formula for the optical rotation 5 can now be written in the form 

s = 16n 2 Nz y v 2 R ah 

3hc Y v lb ~v 2 + Wlab 

where R a b = Im (a | d | b)(b | m | a) is the so-called rotating power. By separating the 
real and imaginary parts of this expression one finds for the actual optical rotation 



and for the circular dichroism 



K 



where the notations from formula Eq. tlMl) are used. The definition of rotating power given 
above corresponds to the dipole-dipole approximation. In the higher order approximation 
the rotating power must be taken in the form R a ^ = Im (a \ d | b)(b \ m | a) + J2 c d( a I d | 
c)(c | Q | d)(d | m | a) . In reality it is very difficult to calculate the matrix elements R a b 
accurately. However, a number of useful approximate formulas have been derived from the 
expressions Eq. (1591) and Eq.( )60|) . For instance, if in some experiment we can see iV peaks in 
the #(A) function and K peaks in the k(X) function, then it is possible to approximate our 
experimental data by using the two following formulas 

where all numerical parameters Aj, Bi, Ci and D{ must be determined by using the experi- 
mental values of 6 and k at different wavelengths A = -. Numerical examples can be found 
in the book by Djerassi jg]. 

Rosenfeld's theory of optical activity allows one to determine the relations between basic 
molecular properties and actual optical rotation and circular dichroism observed in experi- 
ments. Indeed, by applying the known molecular wave functions one can compute the values 
of (a | d | b) and (b \ m | a) which are used in formulas for 6 and k above. By using these 
values we can evaluate the rotating powers R a b- Then we can try to approximate the curves 
9(v) and k{v) obtained in actual experiments. During this step all line widths can be 
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varied as numerical parameters. In practice, this approach works approximately only for 
some simple molecules. For many molecules of interest, e.g., for complex molecules used in 
cancer research, the current accuracy of the numerical determination of the (b | m | a) values 
is not sufficient to make accurate comparisons with experiments. In particular, the signs of 
the (b | m | a) values can be wrong in a number of cases. In addition to this, Rosenfeld's 
theory of optical activity is essentially a semi-classical theory, since all radiation fields in this 
theory are considered classically. The most rigorous analysis of molecular optical activity 
can be performed only on a basis of modern quantum electrodynamics (QED) llj . This 
will be one of our goals in future studies. 



VII. SPECIFIC ROTATION BY CHIRAL ORGANIC MOLECULES 

In previous Sections we have briefly considered some theoretical aspects of the molecular 
optical activity at arbitrary wavelengths. In this Section we discuss a few basic features 
which are known for optically active organic molecules in solutions. In general, if some 
organic molecule has non-zero electric and magnetic moments, then its optical rotation 9, 
Eq.f l59|) . differs from zero. Such a molecule shows a number of phenomena which are usually 
designated as 'optical activity', or briefly, as an optical rotation of plane-polarized light. In 
principle, any molecule which does not coincide with its mirror image can be optically active. 
Moreover, for each optically active molecule one can always find another non-identical form 
of the same molecule which is related to the original molecule through reflection. These two 
forms of one molecule are called enantiomers, specifically D— and L— enantiomers, which 
naturally refer to the right- and left-handed forms, respectively. In many sources the D- 
and L-enantiomers of various molecules are discussed. However, it should be emphasized 
that currently there is no uniform relation between the absolute configuration of complex 
molecules and their ability to be left- and/or right-rotating. Here it is not our intention to 
summarize all basic rules found for numerous organic substances which are optically active. 
Instead, we restrict ourselves to an analysis relevant to observation of optical activity in the 
vacuum ultraviolet region. 

The most interesting cases can be observed in various organic molecules, i.e. in molecules 
which include one or more carbon atoms. Formally, one carbon atom in a molecule which 
is bonded to four different atoms and/or groups of atoms is sufficient for manifestation of 
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optical activity. In general, such a carbon atom is called an asymmetric atom, or a chiral 
center. In many cases the optical activity can be observed in molecules with two, three and 
more asymmetric carbon atoms (for more detail, see, e.g., ^ and references therein). 
A very well known example is the tartaric acid which may exist in the form of D- and L- 
enantiomers and in its meso-form which has no optical activity. Note also that a number of 
organic molecules with no chiral centers show overall optical activity, e.g., allenes, spiranes 
and biphenyls. Such systems are considered as inherently dysymmetric. The active electrons 
in these molecules are delocalized over a chiral nuclear system. 

In this work we restrict ourselves to the consideration of organic molecules with one asym- 
metric carbon atom (chiral center). In general, the observed optical rotation 9, Eq. (l59p . 
produced by one asymmetric carbon atom in a molecule will be small. However, if some 
additional conditions are combined with each other, then the actual optical rotation 9 in- 
creases to moderate, large and very large values. It was shown more than sixty years ago 
that 'close' presence of some special groups of atoms can increase the actual optical rotation 
by a few orders of magnitude. Such special groups are called 'chromophors'. Typical exam- 
ples of chromophors are the — NH 2 , >CO, — CN, — CqR^ groups and some others. It should 
be emphasized that none of these groups is optically active, but each amplifies significantly 
the optical activity of the neithbouring chiral center. 

Currently, there are a few dozens of different atomic groups which are recognized as 
regular basic chromophors and a large number of special groups of atoms which become 
chromophors only at certain wavelengths. In general, any group of atoms which has excessive 
7T— electron density can be considered as a potential chromophor and any experimental study 
of optical activity in organic molecules is reduced to the analysis of various chromophors 
and their influence on one of more asymmetric carbon atom(s). The problem contains many 
complications. For instance, if one of the hydrogen atoms in a benzene ring bonded with 
a chiral carbon center is replaced by the — NH2 group, then one finds the new chromophor 
— C6H4NH2 which has a different influence on this chiral center. In other words, the change 
in optical activity of the chiral center produced by the new chromophor (— C6H4NH2) cannot 
be predicted accurately and uniformly from the analogous information known for the — Cq¥{^ 
and — NH 2 groups. In general, the derealization of 7r— electrons in various organic molecules 
can be used to create a huge number of 'new' chromophors. Moreover, if the same chiral 
center is bonded to both the — CeH 5 and — NH 2 groups, the amplification of its chiral activity 
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will be drastically different from the previous case. In addition to this one also finds that 
substitution of one hydrogen atom in the — CeH 5 chromophor by, e.g., the — CH 3 group 
will also have a noticeable effect on overall optical activity. The following experimental 
analysis must detect (and investigate) a direct relation between the actual optical activity 
and position of the hydrogen atom (in the benzene ring) which was replaced by the — CH 3 
group. 

It is important for the general theory that each chromophore group can be represented by 
a number of poles located in the complex plane of frequencies v (so-called //-plane). In other 
words, any chromophore group has a number of poles in the complex //-plane associated 
with it. If such a pole is close to the real axis, then in experiments one finds large and 
very large values of optical rotation 5 and circular dichroism k, respectively. In real organic 
molecules one finds not one, but a few different chromophore groups. The experimental 
curves for the 5{y) and k,(v) values measured for v > in actual organic molecules are 
the result of interaction between various poles located in a complex frequency plane (//) at 
different frequencies. In general, the interaction between different chromophores produces 
very complicated spectra for the optical rotation 8{v) and circular dichroism k(v). The 
complexity of these spectra rapidly increases as the number of poles per unit frequency 
interval increases. In particular, this is the case for vacuum ultraviolet wavelengths, since 
almost all known chromophors have many absorbtion lines located in VUV region. 



VIII. CONCLUSION 



We have considered the phenomena of optical activity in homogeneous solutions of var- 
ious organic substances. The classical macroscopic theory based on Maxwell equations in 
dielectric (or nonconducting) media is discussed in detail. The Stokes parameters for almost 
monochromatic light are defined rigorously. The relations between the optical rotation and 
circular dichroism are derived from the basic Kramers-Kronig relations. These relations al- 
low one to obtain/evaluate, e.g., the circular dichroism by using the known values of optical 
rotation at the same frequencies. The explicit expression for the tensors of molecular optical 
activity are derived. Our formulas derived for the tensor(s) of molecular optical activity 
can be used to explain a large number of phenomenon currently known in molecular optical 
activity. Note that our formulas can successfully be applied to the case of the Rayleigh (or 
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non-shifted) scattering when U21 = and also to the cases when W21 7^ (shifted or Raman 
light scattering). It is shown that all known lowest order effects of optical activity must be 
described with the use of finite number of tensors (five, three or one tensors). 

We also briefly consider the quantum (or semi-classical) theory of molecular optical activ- 
ity developed by Rosenfeld in [3j . In this theory all molecules are quantum systems, while all 
electromagnetic fields are described by classical Maxwell equations. A possibility to extend 
measurements of optical rotation and circular dichroism into the vacuum ultraviolet region is 
discussed. Currently, this task seems to be extremely difficult, since there are a large number 
of unsolved problems which must be considered before the whole procedure can be usefully 
implemented. Moreover, it is clear that even the requisite experimental technique will have 
many fundamental differences from the technique applied for traditional wavelengths. Nev- 
ertheless, we can expect that measurements of optical rotation and circular dichroism in the 
vacuum ultraviolet region will produce a large volume of very valuable experimental data. 
These measurements will open a new avenue for some important discoveries and improve- 
ments in our current understanding of optical activity of organic molecules. In general, such 
measurements can be performed with the use of high quality VUV-radiation from modern 
electron synchrotrons. 

It should be mentioned in conclusion that many experimental and theoretical methods 
developed earlier to observe/describe the optical activity in solutions of organic substances 
are now a great interest for modern Stellar Astrophysics (see, e.g., Q, [l6] and references 
therein). It was shown that of radiation emitted by many Stellar objects is partially polarized 
In particular, the polarization of radiation emitted by the hot Be-stars is 



(« 1 - 10 %) 

discussed in 17j - [19|. We shall not go into discussion of these and other similar problems 
in this article. Note only that the overall complexity of Astrophysical problems related with 
the polarization of radiation is significantly higher than for 'regular' problems known from 
Organic Chemistry. 
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